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We use ALEPH data on hadronic r decays in order to calculate Euclidean coordinate space cor- 
relation functions in the vector and axial-vector channels. The linear combination V — A receives 
no perturbative contribution and is quantitatively reproduced by the instanton liquid model. In 
the case of V + A the instanton calculation is in good agreement with the data once perturbative 
corrections are included. These corrections clearly show the evolution of a s - We also analyze the 
range of validity of the Operator Product Expansion (OPE). In the V — A channel we find a dimen- 
sion d = 6 contribution which is comparable to the original SVZ estimate, but the instanton model 
provides a different non-singular term of the same magnitude. In the V + A case both the OPE and 
the instanton model predict the same d = 4 power correction induced by the gluon condensate, but 
it is masked by much larger perturbative contributions. We conclude that the range of validity of 
the OPE is limited to x < 0.3 fm, whereas the instanton model describes the data over the entire 
range. 



1. Quantitative understanding of the interface be- 
tween perturbative and non-perturbative effects is the 
central problem in QCD dynamics. Historically QCD 
sum rules based on the Operator Product Expansion 
(OPE) HH) constituted the first serious attempt to de- 
scribe non-perturbative phenomena in QCD. The initial 
application of QCD sum rules to vector and axial vec- 
tor meson lead to very promising results. It was soon 
discovered, however, that not all hadrons are alike @. 
Phenomenology demands that non-perturbative effects 
in scalar and pseudo-scalar channels, both meson and 
glueball, are much bigger than in the vector channels. 
This fact is not reproduced by the OPE but it was real- 
ized that direct instanton effects appear in exactly those 
channels in which non-perturbative effects are large. This 
observation gave rise to the instanton liquid model j^] . 

The available information on hadronic correlation 
functions, both from experimental data, the OPE and 
other exact results was reviewed in ||. Since then, 
the high statistics measurement of hadronic r decays 
t — * i/ r +hadrons by the ALEPH experiment at CERN 
has significantly improved the experimental situation in 
the vector and axial-vector channel. The purpose of this 
paper is to compare these results with theoretical predic- 
tions, both from the OPE and instanton models. In par- 
ticular, we would like to assess the range of applicability 
of the two approaches and put improved constraints on 
the parameters that enter into the calculations. Trans- 
lating the spectral functions measured by the ALEPH 
collaboration into Euclidean coordinate space correlation 
functions will also allow precise comparison of the exper- 
imental data with improved lattice calculations along the 
lines of g. 

2. In the following, we shall consider the vector and 
axial- vector correlation functions Hv(x) = Uu( x )ja(Q)) 
and U A (x) = (fi a (x)jl a (0)). Here, jftx) = qj^q, 

jt°'{ x ) = 97^75 "T 9 are the isotriplet vector and axial- 
vector currents. The correlation functions have the spec- 



tral representation Q 

Hv,a(x) = J ds pv,a{s)D(\/1,x), (I) 

where D(m,x) = mj {^n 2 x)K\(mx) is the Euclidean co- 
ordinate space propagator of a scalar particle with mass 
m. We shall focus on the linear combinations Hy + n,4 
and Hy — Ha- These combinations allow for a clearer sep- 
aration of different non-perturbative effects. The corre- 
sponding spectral functions pv ± PA which are measured 
by the ALEPH collaboration are shown in Fig. |l|. 

In QCD, the vector and axial-vector spectral func- 
tions have to satisfy chiral sum rules. If we assume that 
Pv{s) ~ Pa{s) = above the maximum invariant mass 
s = m 2 for which the spectral functions can be measured, 
then we find that the ALEPH data satisfy all chiral sum 
rules within the experimental uncertainty. However, the 
central values of the sum rules differ significantly from 
the chiral predictions [Q. In general, both py and pa 
are expected to show oscillations of decreasing amplitude 
0. If we set pv{s) — Pa{s) — above an arbitrarily cho- 
sen invariant mass sq this will lead to the appearance of 
spurious dimension d = 2,4 operators in the correlation 
functions at small x. For this reason we have decided to 
use So = 2.5 GeV 2 , which is slightly below the tau mass 
but allows all chiral sum rules to be satisfied. The reader 
should be aware of the fact that we have, in effect, slightly 
moved the data points in the small x region within the 
error bars reported by the ALEPH collaboration. Finally 
we add the pion pole contribution, which is not shown in 
Fig. [j], to the axial vector spectral function. This corre- 
sponds to an extra term n^(ir) = / 2 m 2 _D(m 7r , x). The 
resulting correlation functions Tlv{x) ± Wa{x) are shown 
in Fig. |2| 

3. We begin our analysis with the combination Hy — 
n^j. This combination is sensitive to chiral symmetry 
breaking, while perturbative diagrams, as well as gluonic 
operators cancel. 
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FIG. 1. Spectral functions v(s)±a(s) = 4n 2 (pv(s)+pA(s)) 
extracted by the ALEPH collaboration. 

In Fig. ^ we compare the measured correlation func- 
tions with predictions from the instanton model. These 
predictions are described in great detail in HJDi an< ^ 
the review The main assumption is that the QCD 
vacuum is dominated by strong non-perturbative field 
configurations, instantons. In the simplest model, the 
random instanton liquid (RILM) , the instanton positions 
and color orientations are distributed randomly. The en- 
semble is characterized by two numbers, the instanton- 
anti-instanton density (N/V) = 1 fm -4 and the average 
instanton size p = 1/3 fm. These parameters were fixed 
a long time ago using the requirement that they must 
reproduce the phcnomcnological values of the quark and 
gluon condensates ||. The agreement of the instanton 
prediction with the measured V — A correlation is im- 
pressive and extends all the way from short to large dis- 
tances. At distances x > 1.25 fm both combinations are 
dominated by the pion contribution while at intermediate 
x the p, p' and ai resonances contribute. 

In order to study the validity of the operator product 
expansion we have to study the short distance region in 
more detail. The OPE predicts that the V—A correlation 



function starts with the the following quark-anti-quark 
operators of dimension d = 4 and d = 6 



n v (x)-n A (x) 

2Uo(x) 



-^-rn{qq)x 



+ —a s (x)(qq) 2 \og{x 2 )x 6 + 



(2) 



The value of the dimension d = 4 operator is determined 
by the Gell-Mann-Oakes-Renner (GMOR) relation to be 
(a:/1.66fm) 4 . Using (qq) = -(230 MeV) 3 and the one- 
loop running coupling constant we also estimate the d = 
6 operator as (x/0.66 fm) 6 . This implies that the d = 6 
operator totally dominates over the d = 4 operator. 

This estimate can be checked by considering the value 
of the d = 6 operator as a free parameter and trying to 
extract it from the measured data. Because the d = 4 
operator is so small we use the GMOR value. A simi- 
lar determination of power corrections was already done 
by the ALEPH collaboration using moments sum rules 
0,|2|Jl3|. Nevertheless, fitting the OPE coefficients in 
coordinate space provides important additional insight. 
The results depend on the coordinate range [0, x m ] used 
in the fit, but for x m < 0.3 fm this dependence is weak. 
The result for x m — 0.3 fm, also shown in Figs. 0M is 



U v (x)-U A (x) 
2U (x) 



1.66 fm 



0.66 fm 



(3) 



We find that the size of the dimension d = 6 term agrees 
with the SVZ prediction. However, the range of conver- 
gence of the OPE is only x <^ 0.3 fm. We also note that 
the accuracy of the data in this regime is very poor. 
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FIG. 2. Euclidean coordinate space correlation functions 
Ilv-(2;) ± Ua(x) normalized to free field behavior. The solid 
lines show the correlation functions reconstructed from the 
ALEPH spectral functions and the dotted lines are the cor- 
responding error band. The squares show the result of a ran- 
dom instanton liquid model and the diamonds the OPE fit 
described in the text. 
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We can also check the short distance behavior of the 
correlation function in the instanton liquid. Instantons 
generate the same d — 4 operator in the OPE but the 
nature of the d — 6 operator is different. To leading 
order in the semi-classical expansion there is no radia- 
tively generated a s (qq) 2 \og(x 2 )x 6 operator, but instead 
there is a non-singular (qq) 2 x 6 term. Such terms are 
dropped in standard OPE, but they are present in the 
correlation functions. The numerical value of this term 
is (ir/0.64 fm) 2 , close to the data and the OPE term. 

4. We shall now focus our attention on the V + A 
correlation function. The unique feature of this function 
is that the full correlator is close to the free field result 
for distances as large as 1 fm. This phenomenon was 
referred to as "super-duality" in || . 

The instanton model reproduces this feature of the 
V + A correlator. We also notice that for small x the 
deviation of the correlator in the instanton model from 
free field behavior is small compared to the perturbative 
0(a s /n) correction. This opens the possibility of preci- 
sion studies of the pQCD contribution. But before we do 
so, let us compare the correlation functions to the OPE 
prediction 



ii v (x) + n A (x) 

2U {x) 



1 



n 



1 

384 



<5 2 (G" )V 



2tt 3 

--^-a s (x)(qq) log(a; 2 )izi 6 



(4) 



Note that the perturbative correction is attractive, while 
the power corrections of dimension d — 4 and d = 6 
are repulsive. Direct instantons also induce an 0{x A ) 

correction 1— (-^) a: 4 + . . . fl^-p^l , which is consistent 
with the OPE because in a dilute instanton liquid we 
have (g 2 G 2 ) — 32n 2 (N/V). This term can indeed be seen 
in the instanton calculation and causes the correlator to 
drop below 1 at small x. 

It is possible to extract the value of Aqcd together 
with the power corrections from the data. Because the 
dimension 6 operator is relatively small we have fixed it 
from a joint fit with the V — A correlator. We find 



u v {x) + n A ( x ) 

2U (x) 



= 1 



a s {x) 

7T 



(t 

~ (o.85fm) 



52 fm 



+ 



(5) 



The value of a s (m T ) — 0.35 |l7fl is consistent with other 
determinations pj, but the value of the gluon conden- 
sate term is smaller than the standard SVZ value 0. 
In fact, the data do not show any kind of "dip" and as 
soon as the d — 4 power correction becomes comparable 
to the perturbative correction it is in dramatic disagree- 
ment with the data. Unfortunately, it will be hard to 
improve on this situation even if high precision data that 
cover a larger range of invariant masses in the vector 
channel become available. Within the range of validity 
of the OPE in the V + A channel, x < 0.3 fm, the power 



corrections are always small as compared to perturbative 
corrections. This makes it doubtful that one will ever be 
able to extract the value of the gluon condensate. 
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FIG. 3. Same as in Fig. |[ but with the Yl v (x) - Yl A (x) 
correlator plotted on a logarithmic scale and Ily(i) +Ha(x) 
shown in more detail. 

5. Finally, we address the purely perturbative con- 
tribution to the V + A correlation function, using the 
instanton calculation as a representation of the non- 
perturbative part of the correlation function. This is 
supported by the fact that instantons provide a very ac- 
curate description of the V — A correlator which is free 
of perturbative contributions. The difference between 
the full correlation function and the instanton calcula- 
tion is shown by the squares in Fig. ^. For comparison, 
we also show the full correlation function with only the 
free field behavior subtracted. At short distance, there 
is no difference between the two curves, and both fol- 
low the first order perturbative result a s (x)/n. At larger 
distances Hv(%) + ^a(x) — 2U°(x) starts to drop, but 
the non-perturbatively subtracted correlator continues to 
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grow. This behavior nicely shows the running of a s even 
at moderately large x At x < 0.3 fm the agree- 

ment becomes even better if the two-loop contribution 
is added, but in this case the Landau pole is reached 
earlier. For this reason, the good agreement of the data 
with the one-loop result even for large x > 0.3 fm may be 
somewhat coincidental. The reason one is able to follow 
the pQCD behavior well outside the usual perturbative 
domain is the remarkable degree of cancellation among 
non-perturbative effects. Further high statistics studies 
of this issue using lattice simulations would be very in- 
teresting. 
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well beyond the usual pQCD domain. In the V — A chan- 
nel, on the other hand, there is an opportunity to extract 
the dimension d — 6 (qq) 2 operator from the data. The 
result agrees with the SVZ prediction, but the accuracy 
is limited by the largest invariant mass accessible in r de- 
cays. In addition to that, the instanton model suggests 
the presence of a non-singular d — 6 contribution of the 
same magnitude. Attempts to extract the d = 4 gluon 
condensate operator from the V + A channel fail because 
in the range of validity of the OPE the d — 4 power 
correction remains a small correction to the pQCD con- 
tribution. 

We conclude that the range of validity of the OPE in 
the vector channels is quite small, x < 0.3 fm. This 
means that there is essentially no "window" in which 
both the OPE is accurate and the correlation function 
is dominated by the ground state. Instantons, on the 
other hand, provide a quantitative tool at all distances. 
This is true even though the vector channels, because 
of the smallness of direct instanton effects, are generally 
considered to be the best system to study the OPE. 

Acknowledgments: We thank M. Shifman for valu- 
able comments. The work is partly supported by US 
DOE grant No. DE-FG02-88ER40388. 



FIG. 4. This figure shows an estimate of the perturbative 
part of the V + A correlation function. The solid line is the 
measured correlation function with the free field correlator 
subtracted. The squares show the measured correlator with 
the instanton contribution subtracted and the dashed line is 
the one-loop prediction. 

6. In summary, we have used the high statistics 
ALEPH data on hadronic r decays to calculate Euclidean 
space correlation functions in the vector and axial- vector 
channel. We focussed our discussion on the linear combi- 
nations I1V (x) ± 14^4 (x) . The combination V — A receives 
no contribution from perturbation theory and provides a 
clean probe for chiral symmetry breaking and the quark 
condensate. V + A, on the other hand, allows for a study 
of perturbative QCD and gluonic operators. 

We have compared the two correlation functions with 
the predictions of the random instanton liquid and the 
OPE. The instanton model provides a very accurate de- 
scription of the V — A correlation function for all dis- 
tances. In the V + A channel the instanton model, sup- 
plemented by pQCD corrections with a running coupling 
constant, also provides excellent description of the data 
for x < 1 fm. 

The remarkable degree of cancellation of non- 
perturbative effects in the V + A channel provides a 
unique opportunity to access perturbative corrections 
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